arXiv:1502.02217v3 [math.PR] 5 Nov 2015 
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Existence and Border Cases 
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Abstract 

Bifractional Brownian motion (bfBm) is a centered Gaussian pro¬ 
cess with covariance 

= 2-'^ -\t- , s,tGR. 

We study the existence of bfBm for a given pair of parameters (h, k) 
and encounter some related limiting processes. 
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1 Introduction 

Classical fractional Brownian motion (fBm) G M, with parameter 

H G (0,1], a centered Gaussian process with covariance 

<’(*>.«) :=i(kr+i«r-i«-i>n. (i.d 

is so widely known and used that it needs no further recommendations. The 
remarkable properties of this class of processes are described e.g. in [17, 
Section 7.2] and in [8, Chapter 4]. 
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Houdre and Villa [9] introduced an extension of fBm called bifractional 
Brownian motion (bfBm) as a centered Gaussian process on M with 

covariance 

:= 2-" ((|s|2« + - |t - . (1.2) 

Usual fBm shows up here when K = 1. 

There is one more special case of bfBm directly related to the usual fBm. 
Consider an anti-symmetrized version of fBm, 

V^^\t) :=W^^\t)t>0. 

It is easy to hnd its covariance 

Ry\s, t) = {s + — 1^ — s, f > 0. 

By comparing this formula with (1.2), we see that bfBm 0 < k < 2, 

consists, up to a scaling factor, of the two independent versions of - 

one for positive, another for negative times. ^ 

Houdre and Villa motivate bfBm just by saying that ’’usual fBm seems 
to be a valuable model for small increments [of real processes], but it appears 
to be inadequate for large increments. It is thus very natural to explore the 
existence of processes which keep some of the properties of fBm but also 
enlarge our modelling tool kit”. Marouby [15] conhrmed this deep guess by 
showing how a family of bfBm’s H = |, K G (0,1), naturally appears as a 
limit in Mandelbrot micropulse model (see also [11, Section 14]). On the 
other hand, in [12] bfBm was used for proving new probabilistic inequalities. 
Initially, Houdre and Villa proved the existence of bfBm on M for 

0<H<1, 0<K<1. 

Later on, Bardina and Es-Sebaiy [2] enlarged the zone of existence. Using 
an idea of Lei and Nualart [10], they proved that bfBm exists on M for 

0<H<1,0<K< min {2, |} . 

To the moment when we started this work, it was still unknown whether 
bfBm exists for any other pairs (h, k). We show below in Proposition 3.1 

^We did not find this identification with anti-symmetrized fBm in the literature, al¬ 
though the sister object, a symmetrized fBm appears in [4] under 

the name of sub-fractional Brownian motion in connection to the limiting behavior of 
occupation time of particle systems. 
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that conditions K < 2 and HK < 1 are necessary for the existence of bfBm 
on M+. 

In the zone H>1,0<K<|, which is most difficult for the research, we 
proceed with spectral analysis and trace a new numerical bound between the 
zones of existence and non-existence. We are guided by a guess of D.S. Egorov 
who conjectured that for any hxed H > 1 there exists a positive k(h) < 
such that bfBm exists for (h, k) with any K < k(h) and does not exist for 
any K > k(h). 

2 Existence arguments 

For reader’s convenience, we briefly recall here (and extend) the key argu¬ 
ments from [9] for the case 0 < K < 1, and those from [2, 10, 14] for the case 
1 < K < 2, proving the existence of bfBm. 

2.1 Case Q < K < 1 

The arguments of Houdre and Villa actually have nothing to do with fBm or 
bfBm, as the following statement shows. 

Recall that a Bernstein function is a function / : M_|_ —)■ M_|_ which admits 
the following Levy-Khintchine representation 

poo 

f{\)=a + b\+ / {1 — e~^^)^i{dx), (2.1) 

Jo 

where a,b > 0 are some constants and /i is a measure on (0, cx)) satisfying 
the integrability condition 

Bernstein functions, many examples and their connections to various helds 
of mathematics are discussed in the monograph [16]. Typical examples are 
A —)■ log(l -I- A) and A —)■ A’^ for 0 < K < 1. If 0 < K < 1, the representation 
(2.1) takes the form 

AK ^ K 

r(i-K) 

Proposition 2.1 Let Y(t),t E M. be a centered process with stationary in¬ 
crements and finite second moments 

a{tf := ¥y{tf. 
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Then for any Bernstein function /(■) there exists a process with covariance 

t) ■= f s,teR. 

Remark 2.2 For fBm Y = we have a{t) = |t|“, thus Proposition 2.1 
used with /(A) = A*^ proves the existence of bfBm with 0 < H, K < 1. 

Proof: For /(A) = a + b\ we simply have 

Rf,a{s,t) = b[a{sY + a{tY - a{s - tY] 

= b [EY{sf+ EY{tf-E{Y{t)-Y{ s)Y] 

= 2bEY{s)Y{t) :=2bRY{s,t). 

Therefore, the process Y(t) := \/2/bY(t) solves the problem. 

Let now a = 6 = 0. In view of the formula (2.1), it is sufficient to hud a 
process on M with covariance 

RA{s,t) ■= (l — exp(—a;((T(s)^ + (T(t)^)) — (l — exp(—a;(T(t — s)^)) 

= exp{—xa{t — s)^) — exp(—a;((T(s)^ + o-(t)^)) 

= exp(—X(j(s)^) exp{—xa{ty) [exp(a;((T(s)^ + — a{t — s)^)) — l] 

= exp(—a;(T(s)^) exp(—[exp(2a; (s, t)) — 1] 

= ex 

for any a; > 0. The latter clearly exist along with processes having covariances 

Rri-rA- □ 

2.2 Case 1<K <2 

Following Lei and Nualart [10], consider the real Gaussian process 

pOD 

:= / (1t > 0, (2.3) 

Jo 

where W is an appropriate uncorrelated Gaussian noise. The process 
is well dehned for K G (0,2). By using (2.2), and analogous formula for 

xe(i,2), 

(e-^-l+xAK‘-'rfx, 
it is easy to calculate the covariance 


p(—a;(T(s)^) expi—xaity) -— i?y(s, tY 
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We have 


rEfc!a(s“ + «"-(s + i)'‘), Ke(0,l), 

^x,o(^’ 0 “ ^^ ~ K — 1, 

i §h!l (-*’'-<"+ (*> + «)"). K6(1,2). 

Next, we rescale time by introducing a process 

xi“’"^(t):=xW(|tp«), teM, (2.4) 


which has the covariance 


H 


'^xf{s,t) := cov 


given by 


\ (i«t“ + w“'‘ - (isi"“+i«t")'‘). 

< 2h (In |s| + In |t| - ln(|s| + |t|)) , 


K e (0,1), 
K = 1, 

KG (1,2). 


If HK < 1, consider the usual fBm W^™'\t),t G M, with covariance from 

( 1 - 1 ), 

() = 1 (ki™+itr-- 1 « - * 1 “) 

and, for K G (1, 2), obtain bfBm just by adding up the independent processes 


B(“’‘^)(t) : = 


+ «eE. (2,5) 


For the boundary case K = 2 the integral representation (2.3) does not 
work but we may simply dehne X^\t)^t > 0, as a degenerated random linear 
process with covariance 

:= cov (^Xi^^(s),xj^^(t)) = 2st = {-s^ + (s +1)^) , 

then let again := Xg^^(|tp“), t G M, as in (2.4), and obtain 

B(“’2)(t) ;= 2“^ x^^’^\t) + 2-^/2 t G M, 
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whenever 0 < H < |. 

In another adjacent case K = 1 the bfBm reduces to the classical 

fBm . We hesitate to call it a boundary case because it separates not 
the zones of existence and non-existence but rather two existence zones with 
different properties. 


In the zone 0 < K < 1 the representation (2.5) does not work because the 
signs in the covariance of are opposite to the desired ones. In exchange, 
we have a representation for fBm 


:= t e M, (2.6) 

with independent processes on the right hand side. This is equivalent to 

Since Xq’^^(-) is a smooth process, it becomes obvious that the local prop¬ 
erties of B are the same as those of fBm IT^™^(-), cf. [18, 19]. 

We also see that if 0 < K < 1 and B*^“’’^)(-) exists, then i?(™)(-) exists 
[7], which simply means HK < 1. In Proposition 3.1 we show that HK < 1 is 
necessary for the existence of B^“’^)(-) whatever K is. 

In the following we prefer to work with a modihcation of the processes 
Xq^\ having simpler covariances. For KG (0,1), H > 0, let 


XW(t) : = 

: = 


K 


r(i-K) 

K 

r(i-K) 


vP(«), 


The respective covariances are 

L - 
K 


R7is,t) = = + + 

= |s|'^“ + |i|^“-(|s|^“+|«P“)“, s,f€R. 
Then (2.6) becomes 

V2 (t) = (t) + V¥ (t), te M, 


s, f > 0, 


( 2 . 8 ) 
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(2.9) 


or, in the language of covariances, 

9 _ p(HiK) I r\K p(HiK) 

Finally, notice that an extension to more general processes similar to 
Proposition 2.1 is also possible for the range 1 < K < 2, cf. [14, Theorem 

3.1(i)]. 

3 Necessary conditions 

First of all notice that we must distinguish the existence of bfBm on M+ 
and on M. This is very different from the case of usual fBm where condition 
0 < H < 1 is necessary and sufficient for the existence in both cases. 

Proposition 3.1 If bfBm exists on IR+, then K < 2 and HK < 1. 

Proof: Since the covariance ') has the self-similarity property 

R^’^\cs,ct) = c'^^^R^^’^\s,t), 

we may transform bfBm into a stationary process by letting 

Stationarity of means that its covariance function depends only on the 

arguments’ difference, i.e. 

cov =: R^^’'"\t 2 -Ti), 

where in our case 

= e-“"2-‘^ (^(l + e^H")^-|e"-l|2™) (3.1) 

= (cosh(Hr))^ - 2(2h-i)k |sinh(r/2)|"“‘^. 

By Holder inequality 

0B"\r)\ = \cov{ut^\0),ui’^\T))\ 

= = 1, (3.2) 
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hence, the function |-Rb’'^^(-)I must be bounded and must attain its maximum 
at zero (this is a common property of all stationary processes). 

In our case, when r —)■ +cx), in (3.1) we have the expansions 

2-Kg-HKr (i^g2Hr)K _ ^-KgHKr 

= (1 + Ke-2“"(1 + o(l))) 

= + 2-’^Ke«(’^-2)"(l + o(l)) 

and 

2-Kg-HKr(gr _ i)2hk _ 

= 2-‘'e™^ (l - 2HKe“^(l + o(l))) 

= - HK2i-‘^e(“‘^-^)"(l + o(l)) 

that yields 

^^’■^^(r) = 2-‘^Ke(‘^-2)«"(l + o(l)) + HK2i-‘^e(™-i)"(l + o(l)). 

Therefore, the boundedness of R{-) implies that both conditions K < 2 and 
HK < 1 are necessary for the existence of B (H,K) 

□ 

Another argument for HK < 1 is given in [14, p.626]. 

Proposition 3.2 The following two covariance based necessary conditions 
hold. 

a) If bfBm exists on M, then K < 2 i^- 

b) If bfBm exists on M+, then K < k(h), where 

k(h) := sup |k : sup ((cosh(Hr))’^ - |sinh(r/2)|^™) < l| 

and k(h) < H“^ /or H > 1. 

Remark: We do not have an analytic expression for the function k(-). Some 
values of k(-) are given in Table 5.1 below. 

Proof: 

a) Assume that bfBm exists on M. Since for its covariance we have 
/?b’^^( 1, 1) = 1, —1) = 1, it is true that 

-1 < Rb ""^1, -1) = 2"’^ [2’^ - 22“'^] = 1 - 2(2«-bK^ 

8 



whereas (2 h — 1)k < 1. 

b) Assume that bfBm exists on M_|_. Then the stationary process 
with covariance exists. Then (3.2) yields 

sup ((cosh(Hr))‘^ - (sinh(r/2))'™) < 1 (3.3) 

T>0 

which is equivalent to 

(2 cosh(Hr))’' < (2 sinh(r/2))^™ + 2'', r > 0. (3.4) 

It remains to notice that if (3.4) holds for some value of K, then it holds for 
any smaller positive value of K, since for any a G (0,1] we have 

(2cosh(Hr))“^ < [(2sinh(r/2))'™ + 2«]“ < (2sinh(r/2))'“^ + 2“A 

Finally, if H > 1 and K = H“^, then 

lim ((cosh(Hr))^ - 2(2H-bK (sinh(r/2))'™) = 2'-« > 1. 

Hence, (3.3) fails for K = H“^. Moreover, by continuity arguments, it also 
fails for all K that are sufficiently close to H“^. It follows that k(h) < H“^. 

□ 

The covariance criteria given in this section are quite elementary. They 
take into account only 2-dimensional distributions of the process. In order 
to get sharper results, we need more rehned spectral methods. 

4 Spectral analysis 

4.1 Stationary processes, covariances and spectral den¬ 
sities 

In addition to the self-similar processes let us 

introduce their stationary versions 



[r) 

■= e-“"IF(«)(e"); 

up 

[r) 

:= e-''"/2xW(e") 


[r) 

:= 

^T,k) 

[r) 
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Notice that is one of the well known versions of fractional Ornstein- 
Uhlenbeck process, see e.g. [3, 5]. 

By the definition of we also have 

^ g-K(2Hr)/2^(K) ('g2Hr^ ^ UP{2Ut). (4.1) 

The covariance functions corresponding to these four stationary processes 
are 

^W(r) := cosh(Hr)-22«-i|sinh(r/2)|2«; 

•= 2cosh(Kr/2) — (2cosh(r/2))’^; 

:= = 2cosh(HKr) — (2 cosh(Hr))’^; 

:= (cosh(Hr))’^-2(2«-i)K|sinh(r/2)|2™. 

The basic equality (2.9) transforms into 

2R^^''\t) = (4.2) 


Let us now pass to spectral representations. Recall that by inversion 
formula any covariance function R{-) of a stationary process such that R G 
Li(M) admits a spectral representation 

/ CXD 

R^^f{u)du, r e M, 

•OO 


and the non-negative summable function /(■) is called the spectral density 
of the corresponding process. We denote fx\ fx’^\ fhe spectral 

densities corresponding to the respective covariance functions dehned above. 
Notice immediately that relation Rx’^\^) = -^xH2H'r) yields 



4.2 Spectral criterion for the existence of 

Proposition 4.1 Let KG (0, 1),H > 0. Then bfBm exists on M+ iff 

/?■“>(«.) <2/r'(«). « €R. (4.4) 
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Proof: a) Assume that (4.4) holds. Then by (4.2) 




' —OO 
ROO 


e^’ “ I 2 fw'^^u) - fP''\u)) du 
R^^f{u)du, 


where /(■) is a nonnegative integrable function. It follows that a stationary 
process G M, with covariance R^’^\-) exists, and we obtain bfBm 

by letting 

B(«’«)(t) :=t“’^174“’''^(lnt), t > 0. 


b) Conversely, if a bfBm exists on M+, then a stationary process 
r G M, with covariance exists. Since this function belongs 

to Li{R), there exists a non-negative spectral density / such that 

/ OO 

P^^f{u)du, r G M. 

■CXD 

By (4.2) it follows that for any r G M 

/ OO 

e'"‘(2/r’-/t“’-2'7) («)* 

■OO ^ ' 

= 2R^^'^\t) - R^x''\^) - = 0 . 

Since the kernel of Fourier transform is trivial, we have 

2 /if - /?•“> - 27 = 0, 

Hence, 

2/r’-/r’=27>o. 

□ 


The criterion of Proposition 4.1 becomes meaningful whenever we have 
explicit formulae for the involved spectral densities. They are found below 
in this section. 


11 



4.3 Spectrum associated to the Lei—Nualart process X 

By using the representation (2.3), we obtain 


K 


r(i - K) 7o 


r(i-K) 

g-Kr/2(i _ 


It follows that 
?(k) 


K 


r(i - 

K) 

K 


r(i- 

K) 

K 


r(i- 

K) 

K 



X I 

e-K-/2(i _ - e-^)r-(^+^^dr 

e-Kr/2(^ - - e-^^)e-’^^dv 


-OO 

r*oo 


-OO 

r*oo 


_p-e(’'+-)w-K./2.T _ - 


)g-K^/2(i_g-e 


5((n + r)5((n)cin, 


r(i - K) J_, 

where g{v) := — e“®”). By applying Fourier transform, we obtain 


- r(i-K)n 

K 


K 


e*™ ^(n) ^(n) du 


OO 

OO 


It follows that 


r(i-K) 

/Fw = 




\9{u)?- 


Now we hnd ^(n). By dehnition 

1 /■“ 


?(m) = 


v^i- 
1 




OO 

OO 


J- 

1 /■' 




e-'“"e-K"/2(l -e-"'')dn 
-e-")dr 

1 /“OO 


V27r{—iu — k/ 2) io 

r(-ra-K/2 + l) r(-ra-K/2) 


\X7r{—iu — k/2) 
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We conclude that 


ff(n) = 

Finally, equation (4.3) yields 


K |r(-ra-K/2)p 


r(i-K) 


2% 


= - 


K |r(^-K/2)| 


2 h r(i -k) 


271 


4.4 Spectrum of the fractional Ornstein—Uhlenbeck pro¬ 
cess 

The layout of calculation is very much the same as for the spectrum of Lei- 
Nualart process. Recall that fractional Brownian motion, as a process with 
stationary increments, admits, for H G (0,1), a spectral representation 




U|h+i/2 


where 


Therefore, 


2 r(2H + 1) sin(7rH) 


and 


[/W(r) = e-“"fFW(e") = 




I fp I H-|-l/2 


Rw\'^) = yl 


- 1) (e‘- - 1) 

|^|h+1/2 |^|h+1/2 


= 2 Re 


= 2 i/h Re 


- 1) (e- - 1) 


r*oo 


' —oo 
poo 


^H-hl/2 ^H-hl/2 

-H(^+-) - l) (e*"" - 1) dv 


= ‘^vlRe g{v+ T)g{v)dv, 

J —OO 

where g{v) := e“™ (e*®“ — l). 
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By applying Fourier transform, we obtain 


Rwi'T) = 2|/H^e/ e^'^'^g{u)g{u)du 


' —OO 
POO 


= ‘^yl Re 


(u)r du 


= 2 cos{Tu)\g{u)\^ du 


’ —OO 

roc 


= 2 yl cos(rM)- 




-u] 


du 


' —OO 

poo 


‘^yl 


^iru \9WV + \g{-u)\^ 


It follows that 


fwi.u) = yl (l?(«)P + l?(-«)r) • 


Now we find 'g{u). By definition 

Y roc 


g{u) = 


V^J-. 
1 r 


V^J-. 
1 r 


V^J- 

i 

J- 

t r 


OO 

OO 


V^Z Jo 
t r 


V^Z Jo 
i 


e-^^^g{v)dv 

^-iiu+H)v ('gie- _ 

e-"" (e'"" - 1) dv 
e-"" e'e'"" dv 
r-^ e*’’ dr 

r~^ (cosr + isinr) dr 


TT 


+ 


ITT 


^/^z cos{7rz/2) 2 r(; 2 ) sin(7r2:/2) 

— ^fZj2 cos(7r;2/2) — i sin(7r2:/2) 
zV{z) sin(7r2:) 

zV{z) sin(7r;2) ’ 


where z = Yi + iu. Notice that the integrals that appear after the integra¬ 
tion by parts must be understood as the main value integrals (they are not 
absolutely converging). 
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We infer that 


|2:p|r(;2)P |sin(7r2:)P 
By using trigonometric formulae 

sin(7r2:) = sin(7rH + inu) = sin(7rH) cos(f7rM) + cos(7rH) sin(f7rM) 

= sin(7rH) cosh(7rM) + i cos(7rH) sinh(7rM), 

I sin(7r;2)P = sin^(7rH) cosh^(7rM) + cos^(7rH) sinh^(7rM), 

we may conclude that 


An), . ^ 2 ^/2 _ e-- + e-- _ 

^ |2:P|r(2:)|2 sin^(7rH) cosh^(7rM) + cos^(7rH) sinh^(7rM) 

r(2H + 1) sin(7rH) cosh(7rM) 

2(u^ + H2)|r(H + sin^(7rH) cosh^(7rM) + cos2(7rH) sinh^(7rM) 
An interesting special case is H = 1/2 where, using a formula 6.1.30 from 


we get 


|r(l /2 + m)p = 


A-An) = 


cosh(7rM) ’ 


^ ' 7r(4u^ + 1) ’ 

in accordance with the classical Ornstein-Uhlenbeck covariance function 

'Hr) = 


There is an alternative approach to the computation of the spectral den¬ 
sity H 7 ^ 1/2, due to Barndorff-Nielsen and Perez-Abreu, [3]. Writing 


and using Taylor expansion 


(1 - e-A“). r>0, 




r(2H +1) 

r(A; + l)r(2H -k + l)’ 


one obtains 


^ —hi r A » 

Ri^\T) = — + ^(-1)^+1^ , r 


whereas 


/;?(«) 


k - H 

27r(u^ + (k — h)2) 
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5 Computation 

According to Proposition 3.2, in order to check the existence of bfBm on M_|_, 
one must check for each pair (h, k) with H > 1, 0 < K < H“^, whether for all 

M > 0, 

^ K |r(if-K/2)P 

Jjf l ) 2 h r(l - K) 27r 

< (5.1) 

r(2HK + 1) sin(7rHK) cosh(7rM) 

(u^ + (HK)2)|r(HK + ra)P sin^(7rHK) cosh^(7rM) + cos2(7rHK) sinh^(7rM) 

Our numerical computations show that for every hxed H > 1 there is a 
positive bound k such that condition (5.1) holds for all K < k and does not 
hold for all K > k. 

More precisely, we went through the range H G (1,100] with a step 0.01. 
For every H we went through the range K G (0, H“^) with the same step 0.01 
and checked inequality (5.1) for u G [0, 50] with the step 0.01. Computations 
indicate that larger values of u are by far irrelevant. They also show the 
existence of the boundary k := k(H), separating the existence and the non¬ 
existence zones as said above. Furthermore, we repeated the procedure with 
the smaller step 0.001 in K near the boundary value, in order to get sharper 
values of function k. 

Few words about the boundary case: since the expression for covariance 
function of bfBm is continuous in H and K, it is clear that the property of 
its non-negative definiteness is conserved when we let the parameters go to 
some limits. In other words, the domain of existence of bfBm is closed on the 
plane {H, K). In particular, the bfBm with boundary parameters (FT, K{H)) 
should exist. 

A sample of values of k(H) is given in Table 5.1 along with its upper 
bound k(h). The complete table of values of functions k(H) and k(h) may 
be found in [13]. The bound k(h) is rather sharp, especially for large values 
of H. 

The resulting global zone of existence for bfBm {B*^“’’^^(t),t > 0}, is 
represented at Figure 1. 
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H 

K(H) 

K(H) 

H 

K(H) 

K(H) 

H 

K(H) 

K(H) 

1.01 

0.988 

0.988 

2 

0.422 

0.440 

6 

0.117 

0.123 

1.1 

0.887 

0.894 

2.5 

0.321 

0.338 

7 

0.099 

0.104 

1.2 

0.794 

0.807 

3 

0.260 

0.273 

10 

0.067 

0.070 

1.3 

0.718 

0.734 

3.5 

0.217 

0.228 

20 

0.032 

0.033 

1.5 

0.603 

0.619 

4 

0.185 

0.196 

60 

0.010 

0.010 

1.7 

0.519 

0.533 

5 

0.144 

0.152 

100 

0.006 

0.006 


Table 5.1: Existence boundary k(h) and its covariance upper bound k(h) 



Figure 1: Global zone of existence for bfBm 

6 Some boundary cases 
6.1 A limiting process for K = 0 

Let us consider a limiting behavior of the bfBm covariance function when 
H > 0 is hxed and K —)■ 0. For r > 0 we have 

= K~^ (exp (Kln(cosh(Hr))) - exp (k((2h - 1) In 2 + lnsinh(r/2)))) 

—)■ ln(cosh(Hr)) — (2 h — 1) In 2 — lnsinh(r/2) 

= ln(2 cosh(Hr)) — 2 h ln(2 sinh(r/2)) 

= [ln(2cosh(Hr)) — Hr] + 2h [r/2 — ln(2sinh(r/2))] := Ri{t) + R 2 {t). 
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We want to find the spectrnm corresponding to this limiting covariance^. Let 


/i(w) ^ y e '''"Rj{T)dT, J = l,2, 

denote the corresponding spectral densities. In order to hnd the densities 
fj, we use the classical relation between the differentiation and the Fourier 
transform, 

^iu) 


fjiu) = 


y/^{iu) ’ 


J = l,2. 


Since 


R'lir) = 




gHr g-HT 


Hsgn (Hr) = H(tanh(Hr) — sgn (Hr)), 


we have 


R[{u) = tanh —sgn (u/u). 

Furthermore, since 

tanh —sgn (u) = —i\[T\~j2 (sinh(7ra/2))“^ + u~^^ 

we obtain 

tanh —sgn (u/H) 


/i(w) = 


\/^{i 


lU) 


= — (2nsinh(7ra/2H)) - - := H ^(^{u/u), (6.1) 


where 




^(m) = tt 


1 - 




nu 


2 sinh(7ra/2) 

This is a nice function with hnite limit at zero and quadratic decay at infinity. 


Next, easy calculation shows that for r > 0 


i? 2 ('r) = H (1 — coth(r/2)) 
-2He-^ 


-2h^. 


n=l 


^Notice however a logarithmic explosion of the term R 2 at zero. This means that the 
limiting process is not a usual process defined pointwise but a generalized one. This feature 
may not be repaired by time scaling. 
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Hence, 


whereas 


and 


= -2Hsgn (r) ^ e-’"'"'”, r G M, 

n=l 



Auiu 


n=l 


f2iu) 



2h 

TT 


n=l 


The spectral density / 2 (-) is locally nice but it decays like |m|“^ at inhnity. 
Therefore, it is not integrable and corresponds to a generalized Gaussian 
process. 


By summing up, we obtain the spectral density 


TT OtT 

fiu) = Mu) + f 2 iu) = - (2Msinh(7rM/2H))-' + ^ + — 

TTU^ TT 

n=l 


Presence of hyperbolic functions in the computations suggests that there 
should be some relation of the introduced objects to hyperbolic geometry. 
This is indeed the case. Cohen and Lifshits studied in [6] many random 
helds and processes on the hyperbolic space. In particular, they introduced 
so called quadratic held playing important role in a hyperbolic version of 
spectral representations. As shown in [6, Section 10.1], being restricted on 
a geodesic line of the hyperbolic plane, quadratic held generates a centered 
Gaussian process with stationary increments Z(r), r G M, with the structure 
function 

EZ(r)^ = 21ncosh(r/2), taw G M. 

Notice by the way that the derivative Z'{-) is a stationary process with the 
spectral density 

fz{u) = . -1. y 
2 smh(7rM) 

A similar expression already appeared in (6.1). 

Let us hx TT = 1 and denote a stationary Gaussian process with 

spectral density (6.1). 

Then straightforward calculations show that independent copies of Z{-) 
and B^^’°^(-) are connected by 

Z{2t) + - bS^’°^(0) = V2W{t), 

where W = is a Wiener process. 
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6.2 Case H = 1: integral representation 


Very few is known about white noise integral representations of bfBm (for 
other processes, see e.g. [11, Section 7.3]). We present here one for the 
boundary case H = 1, 0 < K < 1. By using (2.2), we have 


j x-^-^dx 

{2xsty 


_^_/ 

2Kr(i - K) io 

K r 


2><r(i - K) 7o 

oo 

_^y 

2^T{1 - k ) 




nl 


X ^ ^dx. 


It follows that 


K 


2^T{1 - k) 


1/2 oo 

E*' 

n=l 


e-x^\^^-^-K^I^Wn{dx), 


where Wn are independent Gaussian white noises on M+ controlled by Lebesgue 
measure. 


7 A posterior discussion 

One of the referees made some remarks related to a more general context than 
the particular results of this note. We also believe that a wider discussion 
might be interesting to the reader. 

The referee states “Since the introduction of bfBm there is a real problem 
of motivation. Apart from self-similarity, this process enjoys no inherent 
property, and there are lots of processes with the same qualitative behaviors: 
a lot of Bernstein functions and functions a can be used and, with the help 
of of Proposition 2.1, new processes may be introduced”. We basically agree 
with that, although we hnd that self-similarity is quite a strong additional 
feature for picking bfBm from the crowd of possible generalizations of fBm. 
In our opinion, the decisive argument for motivation of the interest would be 
hnding some natural models converging to the studied process. So far, only 
the case H = 7 of bfBm appeared in this setting but, as we noticed (may be 
for the hrst time?), this case is directly related to fBm itself and, therefore, 
may not advocate for the entire family of bfBm’s. 

Back to bfBm, the referee remarks correctly that our existence results are 
not the same for M and M+. This is of course the weakness of our spectral 
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approach (that is focused only on M+) and should be considered as a source 
of open problems. Moreover, considering multivariate case only makes sense 
when the existence of bfBm on the whole M will be understood. Solving these 
problems apparently goes far beyond the means of our note. 

Acknowledgement. We are very grateful to both anonymous referees 
for the careful reading of our note and for their insightful remarks. 
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